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Abstract 

We present the structure function ratios F^/F® and _F 2 C / F® measured in deep 
inelastic muon-nucleus scattering at a nominal incident muon energy of 200 GeV. The 
kinematic range 10 -4 < x < 0.7 and 0.01 < Q 2 < 70 GeV 2 is covered. For values of x 
less than 0.002 both ratios indicate saturation of shadowing at values compatible with 
photoabsorption results. 



For footnotes see next page. 
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Since the discovery by the European Muon Collaboration m 1982 that the structure lunction 
i*2 per nucleon in iron differs significantly from that of a free nucleon [1], nuclear effects on 
structure functions have been extensively studied both experimentally and theoretically [2]. 

The behaviour of the ratio Ff/F? as a function of the Bjorken scaling variable x has 
been measured for many nuclei in the range 0.0035 < x < 0.8. For x < 0.05 the ratio 
is smaller than unity; this depletion is referred to as "shadowing". For 0.05 < x < 0.15 
the ratio is larger than unity by a few percent; this is the so-called "enhancement" or 
"anti-shad owing" region. There is then a fall of the ratio known as the "EMC effect" up 
to x ~ 0.6, followed by a rise due to Fermi smearing thereafter (fig. 1). All these effects 
increase with atomic number A. In the Q 2 range so far explored there is little or no evidence 
of a Q 2 dependence. 

Recent experiments have concentrated on the very small x region. In particular the 
Fermilab E665 experiment found that the ratio _F^ e / F® exhibits saturation of shadowing 
for x smaller than about 0.001 [3]. 

In this paper we present results on the ratios of structure functions F£ / F® for the 
isoscalar nuclei Li and 12 C. The ratio F°/F? was previously measured by the SLAC E139 
experiment [4] and by NMC [5, 7]. The ratio F^/F^ was also measured previously by 
NMC [6, 7]; in that case however it was obtained by combining the ratios F^ 1 / F% [6] and 
F% I F® [5], measured at different muon energies. The data presented here for both Li/D 
and C/D extend to much smaller values of x than those available so far. 

In deep inelastic charged lepton scattering from an unpolarised target the double differ- 
ential cross section per nucleon can be written, in the one-photon exchange approximation, 
as 

d 2 <Ti~, 47R* 2 F 2 (x,Q 2 ) 
dxdQ 2 Q 1 x 

Here F^{x, Q 2 ) is the structure function of the nucleon, x = Q 2 /2Mv is the Bjorken scaling 
variable and y = v/E while E, v and — Q 2 are the incident lepton energy, the energy transfer 
in the laboratory frame and the four momentum squared of the virtual photon, respectively. 
The mass M is the proton mass, the muon mass and a is the electromagnetic coupling 
constant. The function R(x,Q 2 ) is the ratio of the longitudinally to transversely polarised 
virtual photon absorption cross sections. Assuming R to be independent of the nuclear 
mass, the structure function ratios F£ / F® are equal to the corresponding cross section 
ratios. This assumption is supported by data from SLAC [8] and from NMC [9]. 



1 _ „ _ mM. i h _ 2 J^l\ vL ( i+4M 2 x 2 /q 2 \ 

1 V IE + I 1 <3 2 I 2 \ l+R(x,Q 2 ) J 



2 The experiment 

The experiment was performed at the M2 muon beam line of the CERN SPS with the NMC 
spectrometer [10, 11]. The mean incident muon energy was 197 GeV with a ±4 GeV spread 
(rms). 

Three independent triggers were used. Muon scattering angles larger than 12 mrad 
were covered by the standard trigger (Tl). The small-angle trigger (T2) was sensitive to 
muons scattered between 5 and 17 mrad. Both triggers covered the region x < 0.7, the T2 
events having smaller Q 2 and v at a given x value. Events at larger x were accepted by 
Tl only. An additional trigger was installed (T14), sensitive to muon scattering angles as 
small as 1 mrad, extending the kinematic range down to smaller values of a;. In order to 
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Figure 1: Knowledge about nuclear effects on prior to this work (see text). A phe- 
nomenological curve (top) and some recent high precision data sets (below). 
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Figure 2: Kinematic regions covered by the physics triggers. 
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Figure 3: Complementary target setup (see text) 



avoid triggering on divergent beam tracks, only the central part of the beam was used. In 
this way an extension of the kinematic range down to x = 10~ 4 was achieved, more than an 
order of magnitude lower than in earlier NMC measurements. Figure 2 shows the kinematic 
regions covered by the three triggers. 

The data discussed in this paper were collected with the target configuration shown in fig. 
3. This configuration consisted of two complementary target sets, an upstream one with Li 
and D targets and a downstream one with C and D targets. Only one of the two target rows 
was used at a time, simultaneously exposing to the beam one lithium and one deuterium 
target in the upstream set, and one carbon and one deuterium target in the downstream 
one. The row was frequently exchanged with the complementary one where the "heavy" 
target material A (Li or C) and deuterium were interchanged. Geometrical acceptance and 
efficiency corrections thus cancel in the calculation of the ratios, as do the beam fluxes. 
The frequent exchange of the two target sets (every 60 min, i.e. about 150 times during 
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cut was applied. 6 is the muon scattering angle. 
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the whole run) substantially eliminated the effects of possible time dependences. The cross 
section ratio then depends only on the number of events N and the number of nucleons T 
per unit area in the upstream (u) and downstream (d) targets: 



a A 



\ / meas \ a a I a a 

With this setup, two ratios were measured simultaneously: PLi/ap with the upstream 
targets and (Tc/pd with the downstream ones. The targets had a similar number of in- 
teraction lengths for all materials, providing optimal statistical accuracy. In order that 
the acceptance for the heavy targets be the same as for the deuterium targets, the former 
were segmented to cover the same length. The liquid deuterium targets were 1 m long; 
the corresponding lithium ones consisted of three 13 cm long cylinders of 6 Li kept m an 
argon atmosphere in plexiglas containers; the cylinders were equally spaced over aim dis- 
tance. Likewise the carbon targets consisted each of five uniformly spaced 2 cm thick disks 
in a dallite tube also distributed over aim distance. The total target thicknesses were 
17.7 g/cm 2 , 17.4 g/cm 2 and 18.7 g/cm 2 for deuterium, lithium and carbon, respectively. 
The carbon targets were of natural isotopic composition. The deuterium targets contained 
1.5 % H atoms from a contamination by HD molecules, while the 6 Li targets contained 4.5 % 
7 Li. 



3 Analysis 

The data were processed using a chain of programmes [11] which performed pattern recog- 
nition and geometrical reconstruction of the incident and scattered muons and of the sec- 
ondary charged particles. 

Cuts were applied to remove events from regions with rapidly changing acceptance and 
poor vertex resolution (small scattering angle), poor kinematic resolution (small v) and 
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Figure 4: The effects of the hadron cuts, a) The original x distribution dominated by 
the /xe~ elastic peak at small x. b) The ^-distribution after requiring a charged track in 
the spectrometer not identified as an electron or positron by the calorimeter, c) Contour 
plot of the event distribution over pj| /pjj c and x ("he" stands for hadron candidate). The 

isolated island around pj| /pjj c = 1 and x = m e /M = 5.4 • 10~ 4 represents the two-body 
/xe~ scattering events, d) The ^-distribution with the two-body events removed. 

high background from hadronic decays (low scattered muon energy). The region at high y, 
where radiative corrections are large, was also excluded. At this stage, the region x < 0.0008 
had to be discarded because it was dominated by /xe~ elastic scattering events. Table 1 
summarizes the cuts applied for the three triggers; we shall refer to them in the following 
as the "standard cuts". They rely solely on the incoming and scattered muon kinematics. 

The elastic /xe~ scattering events occur at x = m e /M, where m e is the electron mass 
(fig. 4a). In order to remove this background we used a second set of cuts called the "hadron 
cuts". At least one hadron with energy larger than 5 GeV fitted to the muon vertex was 
required. Hadrons, which are the signature of a deep inelastic event, were identified and 
distinguished from fast electrons (signature of an elastic /xe~ event) by means of the ratio 
of energies deposited in the electromagnetic and hadronic parts of the calorimeter H2 [12]. 
However, due to an acceptance gap in H2 near the beam, about 3% of the electrons were 
misidentified as hadrons (fig. 4b). 

An improvement of the /xe~ event rejection was achieved by exploiting the kinematics 
of these events. Elastic /xe~ scattering is a two-body process, where the electron carries the 
full momentum lost by the muon (i.e. the momentum of the virtual photon), while in deep 
inelastic scattering this momentum is shared by several final state particles. We quantify 
this criterion by observing the ratio of the longitudinal momentum of the virtual photon p7, , 
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is visible, with a broad continuum of DIS events elsewhere. This made it possible to apply 
a box cut in the pj| /pj| c , % plane. The x distribution after this cut is shown in fig. 4d. 

The contamination by /xe~ events in the two affected a;-bins was estimated to be less 
than 6% in the final sample, the residue being due to electrons which radiate a real photon. 
The requirement of a fast hadron in the final state removed all coherent (/x— nucleus) and 
quasi-elastic (/x— nucleon) events. Radiative corrections, otherwise very large in this low- 
x region, were thus significantly reduced. Since the additional hadron track fitted to the 
vertex significantly improves the vertex resolution, the 9 cut could be dropped. 

In summary, at least one charged hadron, i.e. a non-muon particle with less than 75% 
of its energy deposited in the electromagnetic module of the calorimeter was required; for 
events in the range 0.0003 < x < 0.0008, the hadron was required to carry less than 2/3 
of the virtual photon momentum. To complete the hadron cuts, a higher minimum v cut 
of 60 GeV was imposed for the C/D targets (see section 4.1). The hadron cuts were only 
applied to T14 data. 

The number of events after applying the standard cuts for Li/D is 0.22 • 10 6 , 0.37 • 10 6 
and 0.12 • 10 6 for Tl, T2, and T14 respectively, giving a total of 0.70 • 10 6 . For C/D we have 
0.14 • 10 6 , 0.26 • 10 6 and 0.07 • 10 6 events in Tl, T2 and T14 respectively, totalling 0.47 • 10 6 . 
After applying the hadron cuts, we are left with 67000 T14 events for Li/D, and 33000 for 
C/D, but in a different kinematic region going down to x = 10" 4 . Most of the suppression 
of the T14 events by the hadron cuts is due to the removal of radiative events. 

In order to correct the measured yields for coherent and quasielastic scattering and 
higher order electromagnetic processes, the cross section ratio was calculated from eq. (2) 
after weighting each event with the radiative correction factor tj = <r l7 / cr meas . These ra- 
diative corrections were computed according to the prescription of Akhundov et al. [13]. 
The procedure corrects for the radiative tails of coherent elastic scattering from nuclei and 
of quasi-elastic scattering from nucleons, as well as for the inelastic radiative tails. As we 
discussed above, the corrections for the coherent and quasi-elastic tails were not necessary 
for the sample of data for which a fast hadron was required in the final state. 

Several input parameters are needed for the evaluation of the radiative correction factors. 
The evaluation of the inelastic higher order processes requires the knowledge of both R and 
i*2 over a large range of x and Q 2 . A fit [14] to the F® results of deep inelastic scattering 
experiments was used in the measured region (Q 2 > 0.5 GeV 2 ). For the extrapolation of F^ 
to lower Q 2 various models were used. That by Badelek and Kwiecihski [15] combining the 
generalized vector meson dominance and partonic approaches was used for the radiative 
correction itself; extrapolations based on the Regge parametrisation by Donnachie and 
Landshoff [16] and one resulting from a phenomenological fit [11] to low energy data in 
the resonance region [17], were used to estimate the systematic uncertainty. The structure 
functions F^ and F^ 1 were obtained by multiplying F® by empirical fits to the present 
cross section ratios. For carbon the SLAC-E139 data [4] for x > 0.4 were included in the 
fit. Since the measured ratio was needed as an input, an iterative procedure was required. 
The function R(x,Q 2 ) was assumed to be nucleus independent and taken from a fit [18]. 

For the quasi-elastic tails, the nucleon form factor parametrisation of Gari and Kriim- 
pelmann [19] was used. The reduction of the elastic nucleon cross section with respect to the 
free nucleon one (quasi-elastic suppression) was evaluated using the results of a calculation 
by Bernabeu [20] for deuterium and carbon, whereas for lithium the result of Moniz [21] 
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Figure 5: The upstream part of the vertex distribution in two bins — a low 9, low resolution 
T14 bin (up) and a high resolution Tl bin (down). The events originate (left to right) from 
a proportional wire chamber POH, the deuterium target, the three slabs of the segmented 
lithium target, and from another proportional chamber, POB. The dashed lines represent 
the contributions from the individual target segments to the overall fit (full line). 

was used. 

Finally, in order to evaluate the coherent radiative tails the knowledge of the nuclear elas- 
tic form factors is necessary. Parametrisations of available data were used for deuterium [22] 
and for the heavy targets [23, 24]. 

For all the above input parameters, alternative parametrisations from the literature were 
used to estimate the associated systematic uncertainties. For a more detailed discussion, 
see ref. [10]. 

With the standard cuts the radiative correction factors tj for the individual targets range 
from 0.33 at small x to 1.17 at large x for carbon. However, the resulting mean corrections 
to the ratios are 15% and 33% at the lowest x bin for Li/D and C/D, respectively, and 
they decrease rapidly with increasing x to less than 2% for x > 0.008. For the hadron cuts 
the corrections to the ratio were everywhere smaller than 1 %. 

The finite resolution of the spectrometer leads to an uncertainty in the position of the 
interaction vertex. Consequently events might be wrongly attributed to a target material 
or might fall outside of the target region. To estimate the number of such events, the vertex 
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lowest and highest 6 bin together with the fitted curves. These were used to determine the 
tails of the vertex distributions. Correction factors which accounted for the wrong target 
association were calculated from the value of 6 of each event. The largest correction to the 
ratio from this source was 1.4%, at the lowest scattering angles. 

The target thicknesses were corrected to account for scattering in the air and argon 
between the lithium and carbon slices, and in the mylar windows of the liquid targets. 
Holes were found in the lithium slabs, mapped and then convoluted with the beam profile 
and thus accounted for. A similar procedure was applied to account for the rounded edges 
of the deuterium target vessels. A correction was also applied to take into account the slight 
non-isoscalarity of the targets. Finally, the effects of the limited spectrometer resolution in x 
and Q 2 (kinematic smearing effects) on the ratio were estimated by a calculation convoluting 
the resolution in the scattered muon energy and scattering angle with the multiple Coulomb 
scattering in the target. The total effect of all these corrections on the ratios was less than 
2%. 



4 Results 

4.1 x dependence 

For both target pairs, the data from the three triggers were analysed separately and the 
results were found to agree in the region of overlap and therefore combined. Figures 6 and 7 
show the results for F^/F® and F^ /F®; the ratios are shown separately for the standard 
and for the hadron cuts. 

We observe good agreement between the ratio determined with the standard cuts and 
with the hadron cuts, where the radiative corrections are very small, thus providing a cross 
check between the two methods. Because of the loss of statistics with the hadron cuts, 
we use the hadron points only for the lowest four x bins which are inaccessible with the 
standard cuts. 

Figures 8 and 9 show the x dependence of the structure function ratios, averaged over 
Q 2 . The systematic errors with the standard cuts are dominated at small x by the uncertain- 
ties in the radiative corrections. These uncertainties were estimated by varying the input 
parameters to the radiative corrections program as discussed above. Other contributions to 
the systematic error include the uncertainty in the vertex smearing and kinematic smearing 
corrections, the latter being dominant at high x. The uncertainty on the measurement of 
the incoming and outgoing muon momenta (Sp/p = 0.2%) was also taken into account. An 
additional systematic error of 0.5% is included over the whole x range as an estimate of the 
uncertainty of the complementary target method used to extract the ratios. The individual 
contributions were added in quadrature to obtain the total systematic error. 

The systematic errors with the hadron cuts have negligible contributions from the ra- 
diative corrections. However, the requirement to detect a hadron in the final state might 
bias the sample due to hadron losses. These may occur due to a reinteraction in the target 
nucleus or to a collision with a nucleus further downstream. 

In both cases a bias to the ratio can result if the effects are different for different target 
materials. The intranuclear reinteraction is more probable in large nuclei. In ref. [25] it is 
shown that the differential multiplicity of hadrons created in copper is suppressed by up to 
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Figure 6: The merged results for F^/F® from all three triggers (full boxes) compared to 
the ratios obtained with the hadron cuts on T14 data (open points). The errors shown are 
statistical only. 
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Figure 7: The merged results for F% / from all three triggers (full boxes) compared to 
the ratios obtained with the hadron cuts on T14 data (open points). The errors shown are 
statistical only. 
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Figure 8: The present result for F^/F® compared to the indirect result obtained by mul- 
tiplying F^/Ff and F^/F? a (90 GeV) with F?/F? and F? a /F? (200 GeV), respec- 
tively [7]. The inner error bars represent the statistical errors, the outer ones the statistical 
and systematic errors added in quadrature. The relative normalisation uncertainty between 
the two results, not included in the figure, is 0.8%. 
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Figure 9: The present result for / ' compared to that of ref. [7]. The inner error bars 
represent the statistical errors, the outer ones the statistical and systematic errors added in 
quadrature. The relative normalisation uncertainty between the two results, not included 
in the figure, is 0.5%. 
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Table 2: The ratio F^/F? and the logarithmic slopes d(F^ j ' F$) / 'd(lnQ 2 ) as a function of 
x from this work. The x value, the Q 2 range and mean Q 2 (in units of GeV 2 ), and the mean 
y are also given. The logarithmic slopes for the first five bins are averaged to increase the 
statistical significance. The normalisation error on the ratios is 0.004 and is not included 
in the systematic uncertainty. 
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0.03 ±0.13 
0.13 ±0.36 
-0.5 ±0.8 
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Table 3: The merged NMC results for the ratio / ' and the logarithmic slopes 
d{F^ I ' F2) I d(lnQ 2 ) as a function of x. Results from the present work are averaged with 
the ones from ref. [7]. The x value, the Q 2 range and mean Q 2 (in units of GeV 2 ), and 
the mean y are also given. The logarithmic slopes for the first five bins are averaged to 
increase the statistical significance. The normalisation error on the ratios is 0.004 and is 
not included in the systematic uncertainty. 
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quark hadronise as it is leaving the nucleus. For this reason, in the hadron method for the 
C/D ratio we apply a minimum v cut of 60 GeV. This brings the methods into agreement. 
The lowest four x bins, which are the only ones in the final data set determined with the 
hadron method, are not affected by this cut and no systematic error was assigned for this. 

As to the collisions with other nuclei, also no correction was made but a detailed study 
was done [10]. Calculating the probability of each detected hadron to have been absorbed 
in flight, we investigated by how much the structure function ratio would have changed if 
we had weighted events inversely proportionally to the survival probability of the hadron. 
Combinations with multiple hadrons in the initial and final states were taken into account 
in three different ways, in particular to account for the hadrons that are missing in the 
final state. As this calculation was not a long chain Monte-Carlo result, and as it gave 
small effects to either side of the uncorrected value, depending on the method used, the 
effects observed were not used to apply a correction, but rather to estimate the systematic 
uncertainty by taking the maximal effect observed on each bin. This is the dominant 
contribution to the systematic error in the hadron method. 

The normalisation uncertainty on the ratios is not included in the errors shown in figs. 8 
and 9. It is due to the uncertainties on the target thicknesses and amounts to 0.4% and 
0.3% for the F^/F? and F?/F? ratios, respectively. 

At small x, the results for the structure function ratios are sensitive to the assumption 
AR = R A -R D = 0. For instance, if AR = ±0.04, then the ratios change by up to ±0.01 at 
small x. 

The data cover the range 0.0001 < x < 0.7 and 0.01 < Q 2 < 70 GeV 2 . Our results show 
the well known pattern of shadowing and enhancement. They extend, however, to smaller 
values of x than previous results, by nearly two orders of magnitude for Li/D and by one 
for C/D (figs. 8 and 9). 

The consistency of the present results with those of ref. [7] was checked using the % 2 and 
the run [26] tests and combining [27] the confidence levels of these two independent tests, 
the first measuring the magnitudes, and the second the signs of the deviations. For the 
two F% I F® ratios, the combined confidence level is 3.5%, whereas for the F^/F® ratios 
it is 7.2% if the first bin from ref. [7] at x = 0.0085, showing a 4.9 standard deviation 
discrepancy, is excluded. 

The results for F^/F® are presented in table 2, and those for F^ / ' F^ combined with 
the results of ref. [7] are given in table 3. 

Figure 10 shows a compilation of small x data on Xe/D, Li/D and C/D, including the 
ratios at x = derived from real photoabsorption experiments [28]. 

The shadowing in lithium is found to saturate at a value of 0.890±0.010(stat.)±0.021(syst.) 
for x < 0.002, determined by averaging the data below x = 0.002. In the same region, the 
shadowing in carbon saturates at 0.851±0.014(stat.)±0.029(syst.). The present data are 
consistent with the findings of the E665 collaboration that the Xe/D ratio [3] also saturates 
for x < 0.002, at a value of 0.72 ± 0.03(stat.)±0.07(syst.). 

These saturation values are in agreement with the results from photoabsorption exper- 
iments. At an incident photon energy of 60 GeV it was found that (cr(~fC — ► X)/a(~fD — ► 
X)) = 0.832 ± 0.010 [28]. Similar data are not available for lithium, but we can estimate 
the ratio by power-law interpolation (as suggested by [29]) between deuterium and carbon. 
A weak energy dependence of shadowing has been observed in Cu [28] and the energy of 
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1(T 4 1(T 3 1(T 2 10 _1 1 

x 

Figure 10: The NMC ratios F^/F? (present work) and F?/F? ( merged data, see text), 
shown with data on _F^ e / F® from the E665 experiment [3]. The points at x = show the 
ratio a(~fA — ► X)/a(fT> — ► X) from photoabsorption experiments at _E 7 = 60 GeV [28]; 
the arrows indicate by how much these points would move under the energy dependence 
assumptions described in the text. Dash-dotted lines indicate the saturation values, the 
shaded areas the errors. The errors on the DIS data points are statistical only; the saturation 
bands include statistical and systematic errors added in quadrature. 

the virtual photon in our data is higher than 60 GeV. The arrows in fig. 10 show how the 
real photon points would be affected by an analogous energy dependence for C and Li. 

4.2 Q 2 dependence 

In order to study the Q 2 dependence of the structure function ratios, the latter were eval- 
uated in bins of x and Q 2 (figs. 11-12). For each x bin straight lines in In Q 2 were fitted 
through the data points in order to extract logarithmic slopes d{F^ I F^) /d(lnQ 2 ). The 
slopes thus obtained are shown in figs. 13 and 14. The slopes from the lowest five x 
bins were averaged to increase the significance of the result. The slopes in each x bin are 
statistically compatible with no Q 2 dependence. 

5 Summary and conclusions 

In summary, we have measured nuclear effects on the structure function F? in lithium and 
carbon, particularly at small x. Compared to earlier experiments, the kinematic range 
is extended down to x = 10~ 4 . The increase of shadowing and antishadowing with the 
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Figure 11: The ratio F^/F? versus Q in x bins. The mean x of each bin is indicated in 
the figure. The diamonds correspond to T14 data analysed with the hadron cuts, and the 
black points represent the results obtained by merging the data from all triggers analysed 
with the standard cuts. The errors shown are statistical. 
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Figure 12: The merged NMC results for the ratio F^ I F^ versus Q 2 in x bins. The mean x 
of each bin is indicated in the figure. The diamonds correspond to T14 data analysed with 
the hadron cuts, and the black points represent the results obtained by merging the data 
from all triggers analysed with the standard cuts together with those of ref. [7]. The errors 
shown are statistical. 
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Figure 13: The logarithmic slopes d^^ 11 / F®) /d(lnQ 2 ) as obtained from linear least squares 
fits to the points of fig. 11. 

atomic number A, as well as the lack of significant Q 2 dependence in the measured range 
is confirmed. In the new kinematic region at very small x, we have found an indication of 
saturation of shadowing at values compatible with the results of photoabsorption experi- 
ments. 
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(this work & ref. [7]) 
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Figure 14: The logarithmic slopes d{F^ / F^) /d(lnQ 2 ) as obtained from linear least squares 
fits to the points of fig. 12. 
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